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We report on a study of the shape of a stiff, charged rod that is subjected to equal and oppo- 
site force couples at its two ends. Unlike a neutral elastic rod, which forms a constant curvature 
configuration under such influences, the charged rod tends to flatten in the interior and accumulate 
the curvature in the end points, to maximally reduce the electrostatic self-repulsion. The effect of 
this nonuniform bending on the effective elasticity and on the statistical conformations of a fluctu- 
ating charged rod is discussed. An alternative definition for the electrostatic persistence length is 
suggested. This new definition is found to be consistent with a corresponding length that can be 
deduced from the end-to-end distribution function of a fiuctuating polyelectrolyte. 

PACS numbers: 82.35.Rs, 87.15.La, 36.20.-r, 82.35. Lr 



I. INTRODUCTION AND SUMMARY 

Given the ubiquity of charged linear structures in biol- 
ogy, and their fundamental importance to essential pro- 
cesses in living systems, the relationship of Coulomb 
interactions to the mechanical properties of polyelec- 
trolytes is a topic of pressing interest. Indeed, there have 
been a number of theoretical P, 0, 0, S IS and ex- 
perimental 0, H, 13 studies of the behavior of charged 
polymeric chains, with an eye to elucidating the vari- 
ous influences that control their equilibrium and statis- 
tical characteristics. In spite of the considerable effort 
expended, there is, as yet, no comprehensive theoretical 
description of the way in which charged chains respond 
to environmental influences. 

The basic theoretical model of a polyelectrolyte chain 
(PE) is simplicity itself. A rod with an intrinsic elasticity 
quantified in terms of a bending modulus carries charges, 
either uniformly distributed along it, or concentrated at 
points along its axis. The energy of this chain consists 
entirely of the elastic energy associated with bending of 
the rod and the electrostatic energy of interaction of the 
charges on the rod. The electrostatic interaction may be 
screened by counterions in solution in the vicinity of the 
rod. This screening is assumed to be Debye-like. 

The characterization of the effects of intrinsic stiffness 
of a neutral inextensible, or worm-like, chain (WLC) in 
terms of a persistence length is by now well-established 
[Toj. This quantity describes the exponential decay of 
correlations in the orientation of the backbone of a WLC. 
It is directly related to the energy stored in a short seg- 



ment of the WLC which has been bent as the result of 
the application of force couples at its end points. An 
extension of t he p ersistence length to the case of a PE, 
due to Odijk |lj utilizes this basic approach to obtain 
that quantity in the case of a stiff, charged rod. In cal- 
culating the energy of the bent segment of the PE, Odijk 
assumed that the shape of this segment is not affected 
by the Coulomb interactions, an assumption shared by 
all other expressions for the effective persistence length 
of a PE of which we are currently aware This 
leads to a remarkably successful analytic expression for 
the persistence length [see Eq. © in the next section]. 
This approximation does not work in all regimes; in pre- 
vious work, we have identified the regimes in which this 
assumption is correct and other regimes in which it fails 
to give rise to accurate answers JJJ . 

The general question of the applicability of the no- 
tion of a persistence length to a PE in the rod-like limit 
was explored in the above mentioned work |l5l Il6| by 
the current authors. This study was conducted in the 
context of a calculation of the thermal distribution of 
end-to-end distances of an ensemble of rod-like PE's. It 
was found that there are regimes in which this distribu- 
tion differs substantially from the corresponding distri- 
bution for neutral WLC's. When this is the case, the 
effective persistence length associated with the end-to- 
end distance distribution (or radial distribution) is not 
consistent with the formula derived in Ref. The ef- 
fective persistence length is obtained in reference [Tsl 
by matching the PE distribution as closely as possible 
to the radial distribution of uncharged worm-like chains 




FIG. 1: The bent rod. The angle 60 between the two ends 
is illustrated. The relationship between the angle Oq and the 
energy, £, stored the bent rod is displayed in Eq. |Q 

In this paper, we undertake a calculation of the effec- 
tive persistence length. Our approach is identical in over- 
all philosophy to the that utilized by Odijk ^llj. A rod 
is subjected to external torques that causes it to bend. 
The energy, £, in the bent rod is calculated and related 
to, ^0, the difference in angles at the two ends, as shown 
in Fig. n The precise relationship between the energy, 
£, in the rod and the angle 9o is 

kBT 2 L' ^ ' 

Here, L is the total length of the rod. In Odijk's ap- 
proach, it is assumed that the bent rod takes the form 
of an arc of a circle. This is true if all of the energy is 
elastic. However, in the case of interest here, a signifi- 
cant portion of the energy may be electrostatic in nature. 
The fundamental improvement over Odijk's calculational 
method is that we determine the actual shape of the bent 
rod, taking into account the effects of screened and un- 
screened electrostatic interactions. We assume that the 
shape taken is controlled by the requirement that, in the 
absence of thermal fluctuations, the rod minimizes its to- 
tal energy. The expression obtained for the rod's shape 
is directly related to the inverse of the Hamiltonian that 
relates the energy of the bent rod to its distortion from a 
straight line. This approach leads to new results for the 
shape of the rod and to alterations in the energy and the 
persistence length as defined in Eq. 

We find that there are, in certain regimes, dramatic 
differences between the shape of a short PE under the 
influence of force couples at its ends and the shape of a 
similarly torqued WLC. In such regimes the curvature of 
a PE is concentrated at its ends, while the WLC distorts 
into a circular arc, as illustrated in Fig. [3 In regimes 
in which the shapes of the PE and WLC coincide, we 
find that the Odijk's formula for the effective persistence 
length is accurate. However, when the shape of the bent 
PE is inconsistent with the assumption of constant cur- 
vature, the electrostatic persistence length derived from 
Eq. Q deviates from Odijk's result. On the other hand, 
this persistence length is compatible with the persistence 
length deduced from the statistics of the conformations 




FIG. 2: The energy-minimizing shape of a charged rod with 
full treatment of electrostatic interactions and end effects 
(plot (a)). Plot (b) corresponds to the energy- minimizing 
shape of a WLC. A segment of a neutral WLC forms an arc 
of a circle while the curvature of a PE segment can be con- 
centrated at its ends. 

of a fluctuating PE described by the radial distribution 
hmction [Hill. 

The rest of the paper is organized as follows. Section 
^introduces the basic elastic model for PE's, followed by 
a general description on how one obtains the equilibrium 
shape of a PE in Sec. IIIII The results for the shape of 
a PE are presented in Sec. IIVI and it is shown in Sec. 
IVlhow a persistence length can be deduced from the PE 
profile. The suggested electrostatic persistence length is 
then compared with Odijk's formula and also the per- 
sistence length that one can extract from the radial dis- 
tribution of the PE in Sec. IVII which is then followed 
by concluding remarks in Sec. IVIII Some details of the 
calculations are relegated to appendices. 

II. MODEL ELASTICITY FOR 
POLYELECTROLYTES 

As noted above, we treat the polyelectrolyte as an in- 
extensible charged rod. The total energy of such a rod is 
the sum of the intrinsic elasticity and the screened elec- 
trostatic interaction energy. We assume that the electro- 
static interaction is screened by counterions that adjust 
more or less instantaneously to changes in the overall 
shape of the polyelectrolyte. Because of the inextensi- 
bility of the PE's under consideration 0,0], we adopt 
Kratky and Porod WLC model to describe the bending 
energy of the chain In this model, polymers are 

represented by a space curve r(s) as a function of the 
arc length parameter s. The total energy of the chain is 
given by 

£ £pO f'' , fdt{s)Y (3 /-^ -.|r(s)-r(«')| 

(2) 

where t is the unit tangent vector. The second term 
in the Eq. ((SJ is the Debye-Hiickel potential, in which 
screening is controlled by the Debye length that is a 
measure of the ionic strength of the solvent. The quan- 
tity P = £b /b"^ is the strength of the electrostatic interac- 
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tion with b the average separation between neighboring 
charges and £b — e^/eksT the Bjerrum length. The 
quantity e is the dielectric constant of the ion-free sol- 
vent. 

The chain is assumed to be sufficiently stiff that ex- 
cluded volume does not play a role. We will consider 
PE's whose length, L, is either comparable with or long 
compared to the intrinsic persistence length, £pQ. In 
both cases, we restrict our consideration to regions in 
which the combination of intrinsic stiffness and repulsive 
strength of the Coulomb interaction keeps the chains in 
their rodlike limit (ijj. We do not take into account the 
fluctuation in the charges localized to the chain and in 
the counterion system that can give rise to attractive in- 
teractions leading to chain collapse [2(tI l2ll | . 

Orienting the tangent vector at one end of the PE, 
t(0), so that it points in the positive z direction, we 
characterize t(s) by two angles of rotation: Ox{s) = 

^rctHUs)/^l-tUs)-tlis)] in the xz plane and 



9y{s) = arctan[iy(s)/y 1 — t^{s) — t^{s)] in the yz plane. 

For polyelectrolytes in the rodlike limit, 9x{s) and 6y{s) 
are small. With the help of the relation r(s) — r(s') = 

du t(u), we are able to expand both the bending en- 
ergy and the screened Coulomb interaction about the 
rodlikc configuration up to quadratic order in d{s) — 
{0x{s),ey{s)) following the Ref. [H [ill. In this case, 
Eq. ^ can be written: 



knT 



'-^ I ds 



2 7o 1 V '^•^ y \ ds 

+ 7: [ [ ds ds'Cis,s') 
2 Jo Jo 

X {ex{s)e4s') + 9yis)ey{s')}, 



(3) 



where the electrostatic kernel is given by 



C{s,s') 



dsids2 



[1 + k{s2 - Si)] 



X {{S2 - Si) [e(s - Si) - e(s - S2)] S{S - S') 

- [e(s - si) - e(s - S2)] [e(s' - si) - e(s' - s^)]} , 

(4) 

with 8(s) the Heaviside step function. 

Equations Q and I^J constitute the expression for the 
energy utilized by Odijk in his calculation of the elec- 
trostatic persistence length. The essence of this calcu- 
lation is to constrain the difference between the orien- 
tation of the tangent vectors at each end of the rod, 
0{L) — 6{0) = 6*0 and then to determine the energy, £ of 
the bent rod with the use of the formulas in Eqs. Q and 
Q. The persistence length, £p, is then given by Eq. 
In his calculation of the total energy of a bent segment 
of PE, Odijk makes the approximation that the segment 
is characterized by a constant curvature jlll |. That is, he 



assumes that the electrostatic interaction does not cause 
the shape of the PE to differ from that of a WLC seg- 
ment. According to this picture, 9x,y{s) are linear func- 
tions of s. One can then obtain an explicit, analytical 
expression for the total energy of a bent polyelectrolyte 
segment, which leads directly to the following prediction 
for the persistence length of such a charged rod through 
Eq. O: 



tOdijk, 



(5) 



where 

^Odijk 



(31^ 
12 

-f- 



1 

kL 
8 



(6) 



In the next section, we sketch a general method of 
calculating the equilibrium shape of elastic charged rods 
under the influence of applied bending forces. This leads 
to a prescription for the calculation of persistence length 
based on the energy of the bent PE. 



III. THE GENERAL STRATEGY FOR FINDING 
THE ENERGY-MINIMIZING SHAPE OF A 
BENT POLYELECTROLYTE 

As noted previously, the principal advance in our ap- 
proach to the bending energy of a PE is that we calculate 
the actual shape of the bent rod. Here, we discuss the 
general approach to this calculation. The energy of the 
rod in Eq. (jJJ can be expressed as the expectation value 
of the energy operator ,7i, i.e. 



1 



L rL 



dsds' 0{s)n{s,s')0{s'). (7) 



ksT 2 Jo JO 
where Ti., the bilinear "energy" operator, is equal to 



n 



^pO^Sis-s')+pC{s,s'). 



(8) 



To obtain the response of a polyelectrolyte to force 
couples at its two ends, we minimize the energy expres- 
sion in Eq. 10 with respect to 9x,y{s), subject to the 
following boundary conditions: 



dxiO) 

9yi0) 

OxiL) 

6y{L) 



0, 
0, 

f^Ox 

0, 



Oy, 



(9) 



where the angles 0ox and 0oy are assumed to be very 
small, in order to keep the chain in the rod- like limit. 
The resulting Euler-Lagrange equations are completely 
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decoupled with respect to both variables Ox{s) and Oy{s); 
therefore, we focus on planar deformations which can be 
characterized entirely in terms of a single angle 9{s). In 
this case, the boundary conditions are simply. 



^(0) 

0{L) 



0, 



(10) 



We enforce these boundary conditions via Lagrange mul- 
tipliers by adding the terms 



{Xod{s)9{s) - XLd{s - L)e{s)} ds, (11) 



to the energy in Eq. {Tj). Then we seek the solution 
that minimizes the energy £ in terms of an eigenfunction 
expansion of the form 



(12) 



n=0 



where ipn{s) are eigenfunctions of the bilinear Euler- 
Lagrange "energy" operator H in Eq. © . 

The solution to the minimization equation, which is 
now in terms of the amplitudes c„, is 



(13) 



where 's are the eigenvalues of Eq. (|SJ| . The values of 
Ej's depend on the three dimensionless parameters, /3L, 
kL, and £po/L. Rotational invariance implies that there 
must be one eigenfunction, ipo{s) — ^/VL, with eigen- 
value eo — 0. This is due to the fact that simply tilting 
the rod, which yields a constant value for 9{s), does not 
change the energy. All other eigenfunctions possess pos- 
itive eigenvalues. Equation (|13|l shows that there is a 
singularity at eo = 0. We can remove the singularity by 
adding a term of the form |0(s)^ to the energy in Eq. 
(UI). The parameter 5 is a small "gap" parameter that 
will be set equal to zero at the end. The quantity g is 
thus added to the denominator of Eq. H13|l , and there is 
no longer a singularity at eg = unless g = 0. 

We adjust the Lagrange multipliers, Aq and Al, using 
the boundary conditions of Eq. IjlUI) . The two equations 
that lead to the results for the A's arc 



^(0) = 4f^^ + E^»(o) 
VL g ^ 



0, 



(14) 

AoV'n(0) + ALV„(L) 



en +5 



(15) 



The above equations reveal that the only possible way 
in which we can obtain finite values for the angles at 
and L is to have Al —Xq as g — > 0. Let us set 



Ao + Al = gA\/X, where A is a constant that is set by 
adjusting boundary conditions. The general solution for 
9(s) in the limit 5 = is, then, 



0is)^A + Xof2Ms)^^^ 

n=l 



ML) 



(16) 



The limiting result of equal and opposite A's makes 
sense if we think of those Lagrange multipliers in terms 
of torques, or force couples, applied at the two ends of 
the PE segment. Given such a picture, we know that 
unless the two torques are equal and opposite, there will 
be an uncontrolled rotation of the segment. It is, in fact, 
possible to set up a calculation of the shape of a segment 
under the influence of such torques. Energy minimization 
yields equations for the angles at the ends of the segment 
that are precisely as given by Eqs. (|14|l and (|15|l . with 
the A's proportional to the torques at each end. 

In the next section we calculate numerically the 
energy-minimizing shape of a charged rod and compare 
it with an arc of circle. This will lead to greater insight 
into the influence of energetics on the classical shape of 
a bent segment of PE. It will also allow us to test the 
fundamental relevance of an energy calculation such as 
the one described above to the persistence length of a 
PE, defined in terms of its conformational statistics. 



IV. NONUNIFORM BENDING OF CHARGED 
ELASTIC RODS 

The expression for 9 given in Eq. (|16|l can also be 
written as 



{s) = A + Xo{K{s,L)-K{s,0)), (17) 



where 



K{s,s')^f:^tlim^^ (18) 



It is readily demonstrated that the quantity K{s, s') 
defined in (|18|l is the inverse of the energy operator in 
Eq. The operator K{s,s') has been calculated with 
the use of a cosine function basis set in Refs. [THl IT^. 
We utilize our previously-obtained results for the inverse 
operator to numerically calculate the quantity 9 in Eq. 
(|16|l . An outline of the construction of the energy in this 
basis set is contained in Appendix IXI 

In the absence of electrostatic interaction, a bent elas- 
tic rod conforms to an arc of a circle, which is described 
by a linear solution for Eq. (|17|l . In Fig. O the 
energy-minimizing shapes of charged rods with kL — 10, 
£po/L = 0.5, and different values of /?L, are compared 
with that of a corresponding neutral chain. As the fig- 
ure clearly illustrates, a 9(s) that depends linearly on arc 
length, s, does not correspond to the minimum energy 
configuration of bent charged elastic rods, which tend to 




FIG. 3: Plots of the energy-minimizing shape for charged 
elastic rods [Eq. JT^] at kL = 10, ipo/L = 0.5 and (3L = 300 
(dashed line), I3L = 1800 (dashed and dotted line), PL = 
7200 (dotted line). The neutral chain forms an arc of a circle 
(represented by a straight line). As PL increases, charged 
chains tends to flatten more in the interior and accumulate 
curvature in the end points. The extreme range of PL is 
intended to clearly indicate the influence of charging on shape. 

flatten in the interior and accumulate curvature in the 
exterior (see Fig. |21 above). This is in part because the 
electrostatic self-repulsion is lower in the end points due 
to the reduction in repelling neighboring charges there. 
Figure|2|also highlights the fact that deviations from con- 
stant curvature become more pronounced as the charging 
strength, f3L increase. This tendency reinforces the no- 
tion that Coulomb repulsion underlies the concentration 
of curvature near the ends of the bent rod. 

We have also investigated the influence of inverse 
screening length, k, and absolute length, L, on the shape 
of a bent PE segment. Figure ^ illustrates the effects of 
a change of k on the arclength dependence of 9. In this 
figure, the bare persistence length £po = 50 nm and the 
charging parameter /3 = 25 nm~^ corresponding to DNA 
are used, and the length of the segment is set equal to 
100 nm. It is apparent that as the screening length in- 
creases the shape of the bent segment deviates more and 
more from an arc of a circle. 

Additionally, we have looked at the consequences on 
PE shape of changes in the length of the segment, keeping 
all other parameters fixed. Figure shows how changing 
the length, L, causes 9{s) to deviate from a straight line. 
As in Fig. ^ we have set £po equal to 50 nm and (5 
equal to 25 nm^^, as corresponding to DNA. The inverse 
screening length k in Fig. El is fixed at 0.1 nm"'^. 

The nonlocal electrostatic interaction seems to favor a 
decomposition of the linear profile into a piece-wise lin- 
ear one, in which the interior takes up a lower curvature 
and the two end-segments in the exterior acquire a higher 
curvature. The relatively sharp changes in the slope that 
result in "shoulders" in the profile take place symmetri- 
cally at positions denoted by Sc and L — Sc, which can 
be defined in terms of the intersection points of the tan- 
gents to the various segments of the profile, as illustrated 
in Fig. El 

The position of the shoulder Sc is a monotonically de- 



FIG. 4: Plot of 9{s)/6o for a charged rod, corresponding to 
various values of the inverse screening length, k. In all plots 
ipo ~ 50 nm, /3 = 25 nm~^ (corresponding to DNA) and L = 
100 nm. The values of the screening parameters are k — 0.5 
nm~^ (thick line), k — 0.1 nm~^ (dashed and dotted line), 
K = 0.05 nm"'^ (dotted line), and k = 0.02 nm~^ (dashed 
line), respectively. 



e(s)/e, 




0.2 0.4 0.6 0.8 1 

FIG. 5: Plot of 6{s)/6o for a charged rod for which the length 
L is allowed to vary. In all curves £po = 50 nm, /3 = 25 nm"'^ 
(corresponding to DNA) and k — 0.1 nm~^. It is evident 
that longer segments behave more like a WLC in that when 
bent they take a shape with a constant curvature. This re- 
flects the influence on the shape of the combination kL, and 
is consistent with the tendencies indicated in Fig. 21 



creasing function of f3L as shown in Fig. [71 The effect 
of (po/L on the Sc is also illustrated in Fig. [71 As seen 
in the figure, the shoulder Sc increases upon increasing 
the intrinsic persistence length of the chain. The de- 
pendence of Sc on screening, on the other hand, appears 
to be more complicated. One generally expects that as 
kL increases, end-effects become less significant, and the 
value of Sc moves toward zero, resulting in a smoothing 
of the curvature along the chain. However, this is true 
only for strong screening. Figure [HI represents the depen- 
dence of Sc on kL in the strong charging regime, where 
we observe that Sc has a relatively weak dependence on 
screening: Sc slowly increases as kL is increased and then 
starts decreasing with further increase in kL. 

One can understand the appearance of the shoulder 
region as an end effect. The nonlocal nature of the elec- 
trostatic self-interaction leads to enhanced repulsion in 
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FIG. 6: The position of the shoulder Sc can be obtained as 
the crossover point of the tangents to the various segments of 
the profile. 
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FIG. 7: The position of the shoulder 
the charging parameter /3L, for kL — 
(crosses), 0.1 (open circles) and 0.5 (filled circles). As ex- 
pected, Sc increases as the rod becomes intrinsically more 
stiflt. 



the interior of the PE as compared to the end-segments. 
This is mainly due to the fact that there are fewer neigh- 
boring pairs at the end-segments to contribute to the 
mutual repulsion. In other words, one might think of a 
crossover length scale at the two ends, below which the 
intrinsic rigidity (that yields a local resistance to bend- 
ing) dominates the energetics of the chain, while beyond 
that length scale (i.e. in the interior of the PE) it is 
the combination of the intrinsic rigidity and the electro- 
static repulsion that controls the energetics. Interest- 
ingly, such a crossover length scale has been introduced 
by Barrat and Joanny in their study of the length-scale 
dependence of the PE rigidity P, 0] . It is important to 
note that the Barrat-Joanny crossover length is defined 
for the crossover in the fluctuations of the angle {6{s)^) 
(that also has a piece- wise linear dependence on s). Since 
the distribution of the angle 9{s) as controlled by the 
Hamiltonian in Eq. Q is Gaussian, both {6{s)'^) and 
the energy-minimizing 9 (Eq. I17|) are linear functions of 
■^(s, s')~^ (Eq. Thus we expect that in general the 

two crossover length scales coincide. Barrat and Joanny 
propose an expression for the crossover length as 



Sc,BJ 



pO 



(19) 



which exhibits the limiting forms of Sc,bj ~ \/ ^po IP for 

«\/^po//3 < 1, and Sc,bj ^ for Ki/^po//? > 1- 

The Barrat-Joanny crossover length shows a qualita- 
tively similar behavior to the shoulder position Sc as de- 
scribed above, except for the slow initial increase in Fig. 
|Sl(for Sc as a function of k). 

There is a way, to reconcile this behavior with the 
Barrat-Joanny picture. This can be achieved by consid- 
ering the fact that in their derivation of the expression 
in Eq. (|19|l above, they have neglected a logarithmic de- 
pendence in the electrostatic nonlocal kernel for technical 
simplicity. While it is not possible to calculate the correct 
crossover length in a compact form as in Eq. H19|) when 
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FIG. 8: The position of the shoulder Sc as a function of the 
screening parameter kL, for £po/L = 0.5 and f3L = 2500. 
The re-entrance behavior is characterized by an initial slow 
increase followed by a relatively faster decay at larger values 
of the screening parameter. 



the logarithmic factor is taken into account, one can ex- 
tract the limiting forms of the augmented Barrat-Joanny 
(aBJ) crossover length as 



Sc.aBJ 



In 



V^po/P for kV^po/I^ < 1. 



for Ky/ipo/P > 1, 



(20) 

which now exhibits an initial increase in qualitative 
agreement with Figs. and |H1 While this picture can 
qualitatively account for the aforementioned behaviors, 
we have not yet been able to achieve a quantitative char- 
acterization of the shoulder position Sc as a function of 
the three dimensionless parameters £po/L, kL, and /3L, 
and in particular compare it with the dependencies as 
suggested by Eq. (|20|l above, due to the insufficiency of 
the numerical data. 

Having obtained the shape of a charged rod numeri- 
cally, we can now follow Odijk and calculate the persis- 
tence length of PE's. 
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FIG. 9: The electrostatic persistence length, £e obtained by 
Eq. 1221 for kL — 10 and Ipo/L = 0.5 at different values of 
/3L, and comparison with Odijk's persistence length as given 
by Eq. 0, and the electrostatic persistence length that can 
be deduced from the radial distribution function of PE's p^ . 

m 



V. DERIVATION OF THE PERSISTENCE 
LENGTH USING ODIJK'S METHOD 

There is a straightforward way to calculate the energy 
of a bent rod, based on the expression for the angle as 
a function of arc length. The quantity 9(s) given in Eq. 
((T7|l is, to within an additive constant, proportional to 
K{s, L) — K{s, 0), and the energy of a bent charged rod 
can be written as 



Details of the calculation leading to Eq. (|21|) are pre- 
sented in the Appendix IbI According to the definition of 
the persistence length in terms of the energy of the bent 
rod, Eq. Q, we have 

" KiO,0)~K{0,L) ~ ^^^^ 

where is the electrostatic persistence length of the 
chain as defined in the Sec. It is important to note 
that the kernel K(^s, s') depends on the parameters PL, 
kL, and ipo/L through the eigenvalues and eigenfunc- 
tions of Eq. . Figure El shows the values of the elec- 
trostatic persistence length, obtained by Eq. (|22|l for 
kL = 10 and £po/L = 0.5 at different values of /3L (tri- 
angles). Odijk's persistence length as given by Eq. © 
is also plotted in the figure for comparison (solid line). 
For small values of /3L, £e coincides with the Odijk per- 
sistence length, fodijk- As PL increases, the deviation of 
£e and ^odijk becomes more significant. The figure also 
displays the value of the electrostatic persistence length 
that one can infer from the distribution of end-to-end 
distances of an ensemble of fluctuating rod-like PE seg- 
ments (open circles) 0| . There will be more on this 
subject in the next section. 

Our general observation from the comparison of with 
^Odijk is that when the two quantities are equal, the bent 



charged rod is close in shape to an arc of a circle. That 
is, 6{s) as a function of s is nearly a straight line as 
in the case of a neutral chain. This indicates that as 
charging increases, end-effects become more important 
and the description of PE's as neutral chains with an ad- 
justed persistence length is inappropriate. It is clear that 
end-effects play a key role in the elasticity of PE's. Such 
effects are also apparent in the statistical conformations 
of the charged rods. 



VI. THE INFLUENCE OF "END-EFFECTS" ON 
THE STATISTICAL CONFORMATION OF PE'S 

The study of the end-to-end radial distribution func- 
tion, G'(r), of a rod-like PE provides an excellent gauge 
of the statistical conformation of polymers. Using the 
expression for energy given in Eq. ^ , we have obtained 
values for the quantity 

G(r)-(^(r-R)), (23) 

where R = r{L) — r(0). The average in Eq. H23() is over 
an ensemble of PE chains. The function G{r) is, then the 
probability that a given chain in the ensemble will have 
an end-to-end distance equal to r 0, . 

With the use of the radial distribution function, we 
have been able to compare the statistical conformations 
of PE's with those of uncharged wormlike chains. 
Figure irni displays the PE end-to-end distribution (solid 
line) along with the WLC distribution (dashed line) in a 
case in which it is not possible to collapse the two distri- 
butions on top of each other. The persistence length of 
the neutral WLC in the figure was adjusted so that the 
location of the maxima of the two distributions are the 
same. The plot of the uncharged WLC is for £p/L = 0.56. 
The distribution is for a PE segment with £pQ/L = 0.01, 
kL — 1, and PL — 360. 

Using these parameters, we also calculated the energy- 
minimizing shape of a PE [Eq. H16() ] as shown in the 
inset of Fig. ^| It is obvious that end effects are not 
negligible in this case and that the response of the PE to 
the bending force is different from that of neutral chains. 
This example indicates a correlation between regimes in 
which the statistical conformations of a PE chain and 
that of a WLC differ and circumstances under which the 
classical, energy-minimizing shape of a PE segment does 
not trace out the arc of a circle. 

There also exist regimes in which the conformational 
statistics of PE chains in the rod-like limit are identi- 
cal to those of WLC's with adjusted persistence lengths 
j2^. For such cases, PE and WLC distributions are in- 
distinguishable to the unaided eye. Figure ITTI illustrates 
an example of this regime. The distribution function of 
the PE with £po/L = 0.0001, nL = 100 and PL = 36000 
completely obscures the distribution function of a WLC 
with the intrinsic persistence length £p/L — 0.876. As 
in the previous example, the energy-minimizing shape of 
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FIG. 10: Comparison of the radial distribution function of 
a PE (solid line) for kL = 1, £po/L = 0.01, and /3L = 360, 
with that of a neutral chain (dashed line). As the inset shows, 
when the two distributions do not match, the equilibrium con- 
figuration of a bent PE is not given by a constant-curvature 
profile. 



G(r) 
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FIG. 11: Comparison of the radial distribution function of a 
PE for kL = 100, £po/L = 0.0001, and /3i = 36000, with 
that of a neutral chain with the intrinsic persistence length 
ip/L = 0.876. The inset shows that when the two distri- 
butions collapse on top of each other, the PE bends with 
constant curvature at equilibrium. 



the PE is shown in the inset. It is clear that the energy- 
minimizing shape of the PE is not distinguishable from 
that of a neutral chain, as the PE also bends with a con- 
stant curvature in this example. 

We have found that whenever there is a virtually per- 
fect collapse of the distribution function of a PE onto that 
of a neutral chain, the persistence length of the neutral 
chain follows Odijk's prediction, in that, ip = £e + ^pOi 
where £p is the effective persistence length of the charged 
chain, and £e — ^odijk uM- noteworthy that in these 
regimes the energy-minimizing shape of PE's is an arc of 
circle in accordance with the approximation utilized by 
Odijk in his derivation of Eq. © . 

Figure ^ compares the values of electrostatic persis- 
tence length, £e obtained by radial distribution functions 
(hollow circles) to Odijk's formula (straight solid line). In 
the figure, the electrostatic persistence length based on 
Eq. H22|l is also plotted (triangles). For small values of 
f3L, £e's obtained through two different noted methods, 
coincide with Odijk's persistence length, ^odijk- As (3L 



increases, the deviation of £e from fodijk becomes more 
significant. However, the electrostatic persistence length 
obtained through the radial distribution function and the 
one found by using the "real" shape of the chain match 
each other quite well. 

As we decrease the quantity kL, the persistence lengths 
obtained by distribution function and energy-minimizing 
shape also start to deviate from each other. This points 
to the fact that replacing a PE chain with a WLC with 
an adjusted persistence length is not well-justified in all 
regimes and that one should use care in the utilization of 
the notion of an electrostatic persistence length. 



VII. CONCLUSIONS 

Our investigation of the equilibrium shape of a bent 
PE has yielded three striking results. The first is the 
fact that the Odijk formula, Eq. 0, for the persistence 
length applies almost perfectly to the case of a fluctuating 
PE when the bent equilibrium segment has a constant 
curvature. This is consistent with one of the fundamental 
assumptions underlying the derivation by Odijk [Tll |. 

The second result is a suggestion for an improved cal- 
culation of the persistence length based on the energy 
of the bent PE. This approach appears to yield results 
in much closer accord with the calculations of the radial 
distribution function of fluctuating PE segments, even 
in regimes in which Eq. © does not work. We find 
that the fundamental tactic of extracting a persistence 
length from the equilibrium energy of a bent PE yields 
excellent predictions for the effective persistence length 
of an ensemble of fluctuating PE's over a very wide range 
of parameters — if, however, one performs a conscientious 
calculation of the actual shape of the bent PE. 

Finally, we are able to characterize the shape of the dis- 
torted PE segment in terms of "shoulder" regions, imme- 
diately adjacent to the end-points of the chain, at which 
the curvature is significantly greater than in the chain's 
interior. It seems highly probable to us that issues of 
PE energetics are intimately connected to the quantita- 
tive features of these shoulder regions. We are not yet 
able to claim complete resolution of the questions asso- 
ciated with the energetics and conformational statistics 
of rod-like PE chains. However, the fact that one can, 
at least in principle, systematically investigate the equi- 
librium properties of a torqued PE segment gives rise to 
the expectation of substantial progress in the characteri- 
zation of the action of the important biomolecules in the 
family of PE's. 
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by the National Science Foundation under Grant No. 
CHE99-88651. One of us (R.G.) would hke to thank the 
group of Prof, de Gennes at College de France for its 
hospitality and support during his visit. 
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APPENDIX A: EXPRESSION OF THE ENERGY 
OF THE BENT PE IN A COSINE BASIS SET 

In this Appendix, we outline the method by which one 
expands the energy of the bent PE in a basis set that 
automatically satisfies the free boundary conditions at 
the end of the rod. We begin by expressing the distortion 
of the rod in terms of the two-dimensional vector a = 



{t.j;,ty). This means that 

(axis), ay (s), I) 



t{.s) 



(Al) 



Using the Fourier representation of the screened Coulomb 
interaction we find 



f^l fc^t. 2 fds' exp{^k.[r(.)-r(.')]} 



(27r)3 fc2 + «2 



(A2) 



The quantity kj_ is the projection of the wave vector k on 
the x-y plane. The quantity kz is the z-component of that 
three-dimensional vector. Next, we use the series expan- 
sion a(s) = V2j2^=o ^'^^ (^T^) ^ appropriate for the 
open-end boundary condition, and assume for simplicity 
that r is, on average, oriented along the z-axis so that 
Ao = ds a(s) = 0. This leads to the following 



representation of the Hamiltonian of the PE rod 

r^-t7E(--)'^" + f^ E A„.A™K„., (A3) 



where 



En 



AL 



X < cos 



— oo 



dkz_ 
2tt 



{n/df 



U21 



— (n — m) 



sin (%^) sin 



^"^^kl 

2 2 



' COS 



-{n + m) 



kz [kz - (n - m)7r/L] 
sin (^) sin ( Y - ^ 



kz[kz — {n + m)-K/L\ 



sin(Y-^)sin(Y-^) 
[kz — m:/L){kz — nnr / L) 



sin(¥^ + f^)sin(Y-^) 
(kz + mr/L){kz — rmr / L) 



(A4) 



are the elements of the electrostatic energy matrix in the 
cosine basis set. It is now sufficient to replace a^{s) by 
Ox{s), and similarly for ay(s). 

A thorough investigation of the energy matrix (Eq. 
(|A3I) ') is given in Refs. [Hill. The requirement that 
the coarse graining length b not exceed the smallest wave- 
lengths appearing in the cosine basis set puts a restriction 
on the size of the matrix energy. If the length of the PE 



is L, this means that the size, N, of the basis set satis- 
fies N < L/b. At no point in our calculations was this 
inequality violated. 

An advantage of the cosine basis set, quite aside from 
automatic satisfaction of the open boundary conditions, 
is that when n and m are large, the matrix elements in 
Eq. ljA3|l are dominated by those for which m = n. This 
refiects the dominance of elastic energy at short wave- 
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lengths. 



APPENDIX B: CALCULATION OF THE 
MINIMUM ENERGY OF A BENT ROD 

In this Appendix, the minimum energy of a charged 
chain that is shghtly deformed about the rodhke config- 
uration with the use of Eq. (|18() is derived. We begin 
with the expression for the angle when the ends of the 
rod have been torqucd: 



0{s) cx K{s,0) - K{s,L). 



(Bl) 



Here the total arclength of the rod is assumed to be L. 
The relationship between the kernel K{s, s') as given by 
Eq. (|18|l and the energy operator 7i(s, s') as given by 
Eq. © is 



ds" H{s, s")K{s", s) = S{s - s'). 



(B2) 



To obtain the proportionality constant in Eq. IjBip . let 

us assume that the angle at s = is —9q/2, while the 
angle at s = L is 0o/2. Then, we have 



^ K{s,Q)-K{s,L) 
2 K{0,L)- K{Q,Q)' 



(B3) 



The above kernel is symmetric, in that K{x,y) = 
K{y,x); furthermore, there is reflection symmetry in the 
looped rod in that if (0,0) = K{L,L). The next step 
is to note that the energy of the rod is the expectation 
value of the energy operator, i.e. 



2 JO 



dsds' 9{s)-H{s,s')0{s'). 



(B4) 



If we plug in the solution ljB3|l for 9{s) in Eq. and 
make use of the relation l|B2|l . we end up with the ex- 
pression in Eq. (|21|) for the energy of the bent rod. 
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